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Abstract. The spectral properties of the adjacency (connectivity) and 
distance matrix for various types of networks: exponential, scale-free 
(Albert-Barabasi) and classical random ones (Erdos-Renyi) are eval- 
uated. The graph spectra for dense graph in the Erdos-Renyi model are 
derived analytically. 

1 Introduction 

Studies of the network structure seem to be essential for better understanding 
of many real- world complex systems |1I2I3] . Among these systems are social [I] , 
economic jS], biological [5] systems or networks sensu stricto [7] like Internet or 
World Wide Web. In the latter case effective algorithms for WWW content search 
are particularly desired. The Google search engine of the network search bases on 
the eigenvector centrality |3l8j which is well known in the social network analysis 
and not different from the Brin and Page algorithm [319] . In this algorithm each 
vertex i of the network is characterized by a positive weight Wi proportional 
to the sum of the weights ■ Wj of all vertexes which point to i, where Wi are 
elements of the i-th eigenvector w of the graph adjacency matrix A 

Aw = Aw. (1) 

The concept of eigenvector centrality allows distinguish between different im- 
portance of the links and thus is much richer than degree or node centrality 
[TU] , The adjacency matrix A of the network with N nodes is square N x N 
large matrix which elements a(i,j) shows number of (directed) links form node 
i to j. For undirected network this matrix is symmetrical. For simple graphs 
(where no multiple edges are possible) this matrix is binary: a(i,j) = 1 when 
nodes i—j are linked together else a(i,j) = 0. The set of eigenvalues (or its den- 
sity PaW) of the adjacency matrix A is called a graph/network spectrum. The 
graph spectrum was examined [11112] for classical random graphs (Erdos-Renyi, 
ER) T3] and investigated numerically for scale-free networks [14\ by Farkas et 
al. 15 16J. The spectra of complex networks were derived exactly for infinite 
random uncorrelated and correlated random tree-like graphs by Dorogovtsev et 
al. [17] . Several other examples of networks properties obtained by studies of 
graph spectra are given in Ref. [15] . 



While many papers refer to eigenvalues of the adjacency matrices A, less 
is known about the spectra of the distance matrices D. In the distance matrix 
D element d(i,j) is the length of the shortest path between nodes i and j. 
On the other hand, whole branch of topological organic chemistry for alkens 
was developed for small graphs which symbolize alkens' structural formula [19] . 
There, not only adjacency A and distance D matrix but also their sum A + D 
spectral properties were investigated. 

The detailed description of the distance matrix construction during the net- 
work growth for the various network types is given in Ref. [20) . Other solutions 
of this problem are also known; an example is the Floyd algorithm [21j . During 
the network growth nodes may be attached to so far existing nodes randomly or 
according to some preferences P. When this preference bases on nodes connectiv- 
ity k, P(k) ex k, the scale-free Albert-Barabasi (AB) [33] networks will appear. 
The pure random attachment (P(k) = const) leads to exponential nodes degree 
distribution. New nodes may bring with itself one (M = 1) or more [M > 2) 
edges which serve as links to pre-existing graph. For M = 1 the tree-like struc- 
ture appears, while for M > 1 the cyclic path are available. Let us recall that 
degree distributions ir(k) are 7r(fc) oc & -7 , 7r(fc) tx exp(— k) and Poisson's one for 
AB, exponential and ER networks, respectively [11213] , 

Here we study numerically^] the graph spectra pa (A) for growing networks 
with exponential degree distribution for M = 1 and M = 2. We check the 
eigenvalue density pc(A) of the distance matrix D for AB, exponential and ER 
graphs. In literature known to us these spectra was never checked before. 

The graph spectrum pa(A) for dense graph in the ER model is derived ana- 
lytically in Sec. 12. li as well. Here we profit much from Ref. [2"2"] . 

2 Results and Discussion 

Here we show densities of eigenvalues p(A) for matrices A and D for various 
kinds of networks. Results are averaged over N IUD = 10 2 realizations of networks 
of N = 10 3 nodes. 

2.1 Spectral Properties of Adjacency Matrix 

For the adjacency matrix of ER, the density of eigenvalues consist two separated 
parts: the Wigner-semicircle centered over A — and with radius approximately 
equal to 2y/Np(l — p), and the single Frobenius-Perron principal eigenvalue 
near Np [11112123] (see Fig. [Ha)). 

The detailed study of graph spectrum for AB graphs may be found in Ref. 
|15|16] by Farkas et al. There, the deviation for semicircular law was observed and 
Pa (A) has triangle-like shape with power law decay [15]. A very similar situation 
occurs for the exponential networks, but pa (A) at the top of the "triangle" is 
now more rounded. The separated eigenvalues are not observed for this kind of 
networks (see Fig. QJb-c)). 



1 with LAPACK procedure http://www.netlib.Org/lapack/double/dsyev.f 
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Fig. 1. Density of eigenvalues /»a(A) for adjacency matrices A for (a) ER, (b) 
AB and (c) exponential networks with N = 10 3 . The results are averaged over 
iVrun = 100 simulations and binned (AX = 0.1). The isolated peaks in Fig. [Ha) 
correspond to the principal eigenvalue. 



Let us discuss the spectrum of eigenvalues of adjacency matrices of dense 
graphs in the ER model [53] . The diagonal elements of these matrices are equal 
zero a(i,i) = while the off-diagonal elements a(i,j) assume the value 1 with 
the probability p or with the probability 1 — p. The elements a(i,j) above 
the diagonal are independent identically distributed random numbers with the 
probability distribution P(a(i,j)) — (1 — p)S(a(i,j)) +p5(l — a(i,j)). This prob- 
ability distribution of a(i,j) = x has the mean value: Xq — (x) — p and the 
variance a 2 = (x 2 ) — (x) 2 — p(l —p). The universality tells us that the spectrum 
of random matrices does not depend on the details of the probability distribu- 
tion but only on its mean value and variance!! the eigenvalue spectrum in the 
limit N — > oo is identical for different distributions as long as they have the 
same mean and variance. In particular one can take a Gaussian distribution: 
l/V^na 2 exp [—(a; — xo) 2 /2er 2 ] . Thus one can expect that the spectrum of adja- 
cency matrices of ER graphs can be approximated for large N by the spectrum of 
matrices with continuous random variables which have the following probability 
distribution: 



n 
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For the diagonal elements the distribution has the mean equal zero to reflect 
the fact that the corresponding adjacency matrix elements a(i,i) = 0. The last 
expression can be written in a compact form: 
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where 



DA = (V2tt) 
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is the standard measure in the set of symmetric matrices. The matrix B is 
obtained from A by a shift B = A pC where C has the form: 
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The spectrum of the matrix B is given by the Wigner semi-circle law 
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2 if the variance is finite 



It has a support [— 2^/Na, 2yNa], where a = y/p(l — p) as we calculate above. 
We want to determine the spectrum of A. It is a sum A = B + pC of ma- 
trix B for which we already know the spectrum ([5]) and of matrix pC whose 
spectrum consists of an (N — l)-degenerated eigenvalue —p and one eigenvalue 
p(N — 1). The low eigenvalue — p mixes with the eigenvalues of B leaving the 
bulk of the distribution ([5]) intact while the eigenvalue p(N — 1) adds to the dis- 
tribution a well separated peak in the position p(N — 1) far beyond the support 

2y/Np(l — ~pj, 2^Np(l — p) of the main part of the distribution: 

PA (X)*p B (\) + ^5(\-p(N-l)). (6) 

The considerations hold as long as p is finite. For sparse graphs p ~ 1/N — ► 
one sees modifications to the presented picture [22]. 

We note, that the matrix C is both the adjacency and distance matrix for a 
complete graph. Thus two very sharp peaks at A = — 1 and A = N — 1 constitute 
a complete graph spectrum 

PcW = — ^— <5(A + 1) + 1*(A - (N - 1)). 



2.2 Spectral Properties of Distance Matrix 

Spectra of the distance matrix p_d(A) of growing networks for trees (M = 1) and 
other graphs (M > 1) are quantitatively different. For trees the part of spectrum 
for A > is wide and flat. Moreover, the positive and negative eigenvalues are 
well separated by a wide gap [see Fig.[^b-c)] which increases with networks size 
N as presented in Fig. [31 On the other hand, we do not observe any finite size 
effect for negative part of the spectrum. 

The density of negative eigenvalues of D [see Fig. [2] is very similar for consid- 
ered networks. The positive value part of the spectrum for growing networks does 
not depend on growth rule and it is roughly the same for AB and exponential 
networks. 

For complete graph D = A = C and graph spectra consist two sharp peaks 
as mentioned earlier. 



3 Summary 

In this paper the spectral properties of the adjacency A and distance D matrices 
were investigated for various networks. 

For ER and AB networks the well known densities of eigenvalues pa (A) were 
reproduced. For the growing networks with attachment kernel P(k) oc const(fc) 
the graph spectra are similar to the AB networks except of the spectra center. 
For the complete graph two well separated peaks constitute the graph spectrum. 

The spectra of distance matrix D differ quantitatively for trees and other 
graphs. In case of trees (M — 1) the density of positive eigenvalues is very well 
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Fig. 2. Density of eigenvalues pd for distance matrices D for (a) ER, (b) AB and 
(c) exponential networks with N = 10 3 . The results are averaged over iV run = 100 
simulations and binned (Z\A = 0.1). The graphs are horizontally shifted by M 
or p for better view. 
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Fig. 3. Density of eigenvalues pn(A) for distance matrices D for AB trees with 
various network size N. 



separated from the part of the spectrum for A < and extremely flat. Thus the 
specific shape of the distance matrix spectrum may be a signature of absence of 
loops and cyclic paths in the network. 
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